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We present some new three-loop results in lattice gauge theories, for the Free Energy and for the Topological 
Susceptibility. These results are an outcome of a scheme which we are developing (using a symbolic manipulation 
language), for the analytic computation of renormalization functions on the lattice. 
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1. INTRODUCTION 

The computation of vacuum expectation values 
of composite operators is one of the most impor- 
tant tasks in lattice gauge theories. Examples of 
such expectation values in QCD are the Conden- 
sates which are essential tools for the SVZ sum 
rules or the Topological Susceptibility which 
is necessary to solve the so-called C/a(1) problem 
and to understand the 77' mass . In this talk 
we present some results related to the computa- 
tion on the lattice of the Gluon Condensate 



^2 = (0| : ■^Gl^iyd^'^ii^ ■■ |0}, 
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and the Topological Susceptibility x 

X = J d^x{0\T{Q{x)Qm0), (2) 

where Q{x) is the Topological Charge density 
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n(x) = — fl^'^P'^Qb Qb 



(3) 



In these equations, g is the coupling constant of 
the QCD lagrangian and G^^^ the strength tensor 
for the gluon fields. 

To evaluate one of these quantities on the lat- 
tice, e. g. (0|^|0), one defines first a lattice ver- 
sion of A in such a way that A^ a'^A + 
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0(a'^+'^) where d is the mass dimension of the 
operator A and a is the lattice spacing. For the 
Gluon Condensate and Topological Susceptibility 
respectively, our choices for these lattice versions 
are 
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where (x) is a lattice version of the Topological 
Charge density Q 
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In these equations, 11^,^ is the usual plaquette in 
the — V plane. 

Once the lattice version of the operator is de- 
fined, one can perform a Monte Carlo simula- 
tion. The Monte Carlo data will give the physi- 
cal continuum value of (0|^|0) modified by some 
a— dependent renormalizations. These expres- 
sions for the Gluon Condensate and Topological 
Susceptibility can be written in the following way 
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ri>l ^ 
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n>2 ^ n>3 * 
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where N is the number of colors and as usual P = 
2N/g^. The lattice spacing a and f3 are related 
by the renormalization group equation 
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where is the renormalization group invariant 
mass parameter of QCD and tq and ri the first 
two coefficients of the f3 function, /3{g) = —rog^ — 
ng^ - ... 

The last terms in these expressions (those pro- 
portional to c„ and dn) are the perturbative tails 
and represent mixings with the unity operator. 
The second term in Eq.(7) is a mixing with the 
Gluon Condensate. Finally, Zq and Zq are mul- 
tiplicative finite renormalizations which relate the 
lattice and the continuum definitions of the re- 
spective operators. All of these coefficients can 
be computed in perturbation theory and their 
knowledge is essential to extract the physical val- 
ues G2 and X from the Monte Carlo data and 
Eq.(6-7). The values of these coefficients also de- 
pend on the lattice versions used for the opera- 
tors and the lattice action chosen. Throughout 
this work we have used the lattice versions shown 
in Eq.(4) and the Wilson action. 

2. RENORMALIZATION CONSTANTS 

The first coefficients in the renormalization 
terms of Eq.(6-7) can be calculated with rather 
small effort because they involve very few Feyn- 
man diagrams. The result for SU{3) is |]6|-p^ 

Zg = 1 + 0(l//32), 
Zq = l-5.45//3-f 0(1//32) 
b2 = 6.32 X 10^3 ^3 ^ 3 58 ^ IQ-^ 
ci = 2.0 C2 = 1.22. 
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However, the next coefficients involve many more 
diagrams. For instance, Ci and C2 are calculated 
with three diagrams, but the next order, C3, needs 
the computation of 30 three-loop Feynman dia- 
grams! To evaluate these Feynman diagrams we 
have developed an algebraic computer program. 
Major tasks in this algorithm are: 
i) Computing n-point vertices. 
a) Producing a list of relevant diagrams, with the 
corresponding weights. 



in) Performing the contractions for each diagram, 
using up all existing symmetries to produce a 
compact result. 

iv) Extracting the analytic dependence of each 
diagram on its external momenta p, in the limit 
ap 0. 

v) Producing the optimized code for the numeri- 
cal calculation of the loop integrals. 

Regarding the first step, some difficulties inher- 
ent to the lattice are: The existence of vertices 
with an arbitrary number of gluons, a plethora 
of "tensor" structures (due to lack of rotational 
invariance) and a great proliferation in the size 
of vertices (a 6-point vertex in its most compact 
form may typically require some dozens of output 
pages). This in turn necessitates, in the third 
step, simplifying all intermediate expressions as 
much as possible, by devising algorithms which 
use up the symmetries of the diagram under ex- 
change of external legs, under allowed redefini- 
tions of momenta and under permutation of the 
(numerous) dummy indices. 

The fourth step is necessary for computing mul- 
tiplicative renormalizations. The evaluation of 
these form factors is in progress. 

Now the final expression for the diagram can 
be integrated. With the code produced by the al- 
gorithm, we can calculate the numerical value of 
the diagram for rather small lattices and then ex- 
trapolate for larger lattices. This extrapolation is 
performed by assuming the following dependence 
of the diagram on the lattice size L 



diagram = A 



B 
I" 
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The explicit values for A, B as well as that of n 
(which may vary for different diagrams) are de- 
termined by the extrapolation. The result is then 
confronted with the results for an infinite lattice, 
which we also compute. We will explain this com- 
putation for an infinite lattice with an example. 
Let us consider the following integral 



/ = 



"11 1 d^p d\ 



(11) 



which is the value of a two-loop diagram for an 
infinite lattice. In this expression, jp' — 2^^(1 — 
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Table 1 

Exact and fitted values for C3 and 





SU{2) 




SU{?,) 




Exact Fitted 


Exact 


Fitted 


C3 


0.16 0.18 


3.12 


3.0 




7.0 X 10-5 1.8(6) X 10--* 


8.4 X 10-'' 


1.6(1.0) X 10--'' 



cosp^). Using the Schwinger representation for 
the propagators, 

-=/ e--^da, (12) 
P Jo 

we get 

OC CXD 00 

/ = y dai J da2 J da^ $''(ai, Q!2, 013) (13) 

000 
where $ is 
$(ai,a2,a3) = 

J ~7T 27r J— 7r 27T 

exp(2ai cosp + 2a2 cos q + 2q!3 cos(p + qj) = (14) 

-2(ai+Q2+Q3) f + '^dp^ 



g2a2Cosp 1^(2^^ 



2aia3 cosp). 



In Eq.(14) Iq is the modified Bessel function. 
Therefore we have reduced the number of integra- 
tion variables from 8 in Eq.(ll) to 4 in Eq.(13). 
For three- loop integrals with 12 integration vari- 
ables, this procedure reduces this number to 6 
integrations (for a few cases it only reduces to 
8 integrations). Now, a numerical integration is 
feasible. 

The extrapolation done with Eq.(lO) is enough 
to obtain the correct values of the coefficients 
with three digits. Further details of the algorithm 
will be published elsewhere. 

3. RESULTS 

Summing up the individual contributions of all 
diagrams we obtain for the gauge group SU (N) 
and an infinite lattice [pl| 

73 83 

d4 - A^^(iV'-l)(1.73A^2_^0.83+-^^)10-^(15) 



C3 = N^{N^-l){7.28N^-27.35+^^)10-^.(16) 

The exact values for the coefficients obtained with 
our scheme can be compared with the values ex- 
tracted by a best fit of all Monte Carlo data with 
Eqs.(6-9). Both values are shown in table 1 for 
the coefficients d^ and C3 and for an infinite lattice 

The agreement between fitted and exact results 
is manifest for C3. The data for the Topological 
Susceptibility had a rather low statistics, there- 
fore and within the errors the numbers shown in 
table 1 for d^ are in acceptable agreement. 

Finally, we can add the exact values of C3 and 
d4 to those of Eq.(9) and perform again the best 
fits. The values extracted for G2 and x a-i'e in 
complete agreement with those reported in refs. 

MSB 

4. ACKNOWLEDGMENTS 

We wish to thank Adriano Di Giacomo for use- 
ful conversations and the spanish-italian "Accion 
Integrada/Azione Integrata" number A17 for 
financial support. B. A. also acknowledges a 
Spanish CICYT contract. 

REFERENCES 

1. M. A. Shifman, A. I. Vainshtein and V. I. 
Zakharov, Nucl. Phys. B147 (1979) 385, 448, 
519. 

2. G. t'Hooft, Phys. Rev. Lett. 37 (1976) 8; 
Phys. Rev. D14 (1976) 3432. 

3. E. Witten, Nucl. Phys. B156 (1979) 269. 

4. G. Veneziano, Nucl. Phys. B159 (1979) 213. 

5. P. Di Vecchia, K. Fabricius, G. C. Rossi and 



^Actually, for the Topological Susceptibility we calculate 
the unrenormalized value Xu because the physip 
better obtained by using the cooled data of ref. [ 



4 



G. Veneziano, Nucl. Phys. B192 (1981) 392; 
Phys. Lett. BIOS (1982) 323. 

6. M. Canipostrini, G. Curci, A. Di Giacomo 
and G. Paffuti, Z. Phys. C32 (1986) 377. 

7. M. Campostrini, A. Di Giacomo and Y. 
Giindiic, Phys. Lett. B225 (1989) 393. 

8. M. Canipostrini, A. Di Giacomo, H. 
Panagopoulos and E. Vicari, Nucl. Phys. 
B329 (1990) 683. 

9. B. Alles and M. Giannetti, Phys. Rev. D44 
(1991) 513. 

10. M. Campostrini, A. Di Giacomo and H. 
Panagopoulos, Phys. Lett. B212 (1988) 206. 

11. B. Alles, M. Campostrini, A. Feo and 

H. Panagopoulos, Lattice Perturbation The- 
ory Done Symbollically: A Three-Loop Re- 
sult for the Topological Susceptibility, Pisa 
preprint IFUP-TH-31/92; The Three-Loop 
Lattice Free Energy, Pisa preprint IFUP-TH- 
32/92. 

12. M. Campostrini, A. Di Giacomo, Y. Giindiig, 
M. P. Lombardo, H. Panagopoulos and R. 
Tripiccione, Phys. Lett. B252 (1990) 436. 



